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^ : Abstract 
O 

' We present analytical approximations for the real Kelvin function berx and 



the imaginary Kelvin function beix, using the two-point quasifractional ap- 
' proximation procedure. We have applied these approximations to the cal- 

culation of the current distribution within a cylindrical conductor. Our ap- 
proximations are simple and accurate. The infinite number of roots is also 
obtained with the approximation and the precision increases with the value 



o 



o 



^ ■ of the root. Our results could find useful applications in problems where 

analytical approximations of the Kelvin functions are needed. 



I. INTRODUCTION 

Within the scope of special functions, Kelvin functions appear early in the literature 
. They are derived from the Bessel functions of complex argument and many formulas 
can be found. In addition to their importance in Mathematics they have, also, a broad 
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spectrum of application in Physics. For instance, in the analysis of the current distribution 
in cylindrical conductors due to wave propagation |Q. A more interesting application is 
the problem of determining the equivalent impedance of a cylindrical conductor which is 
given in terms of ber x and bei x. Kelvin functions also arose in the problem of the wave- 
impedance of cylindrically layered conductors for application in nondestructive testing 0. 
The temperature distribution in cylindrical conductors due to an alternative current also 
involves berx and beix P]. Kelvin functions also appear in other fields (i.e., fluid mechanics) 
but in this paper we will concentrate on the first electromagnetic application mentioned 
above. 

One of the problems we have found in the applications of Kelvin functions is that their 
behaviour is rather complicated. Kelvin functions are oscillatory and the amplitude of the 
oscillations increases rapidly. In spite of the already existing tables, still the functions are 
part of integrands where the tables are difficult to use. The programs to compute these 
functions are also too elaborated. 

It seems convenient to have a good approximation to the functions, but the conventional 
approximation methods do not give good results because of the pathological behaviour of 
the functions. However, in this paper it is shown how those difficulties can be surmounted 
by using the recently published method of two-point quasifractional approximations [§. 
This has allowed us to find analytic approximations to Kelvin functions with good accuracy, 
and furthermore, the same analytic approximation is valid for the whole range of positive 
values of the variable. 

Power series and asymptotic expansions are required in order to obtain these approxima- 
tions, thus we start analysing those expansions to proceed further away. This will be done 
in Sec. 2, where the form of the approximation will be also derived. Later, the material of 
this paper will be arranged as follows: the calculation of the parameters for the two-point 
quasifractional approximation of ber x is carried out in Sec. 3. Sec. 4 devotes its attention 
to the analysis of the results of the approximation. The same procedure applied to ber x in 
previous sections is employed in the calculation of bei x. That calculation and the discussion 
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of tlie results obtained for beix are performed in Sec. 5. In Sec. 6 we show an immediate 
application to engineering. Finally, Sec. 7 aims its attention to discussion and conclusions 
of the paper. 



II. TWO POINT QUASIFRACTIONAL APPROXIMATION FORM FOR THE 

REAL KELVIN FUNCTION 

As it is well known, the real Kelvin function is defined as 

her X = ^{Jo{iVix)) (1) 

where Jo{z) is the zeroth order Bessel function. 
The power series is given by 

1-^ + ^-- (2) 

The asymptotic expansion of ber x has been also given in several references. Our interest 
here is mainly in the leading term, which can be written as 

e^/^(^cos(^) + ^sin(^)) 
berx ~ 1= —. (3) 

\'27TX 

Once the potential series (Q) and the asymptotic expansion (|^) have been obtained, it is 
possible to determine the form of the quasifractional approximations of the function herx, 
which has an essential singularity at infinity. That kind of singularity is characteristic of 
all the hypergeometric confluent functions. The asymptotic expansion (j^) picks up this 
singularity and shows it up through a branch point at infinity together with essential sin- 
gularities at infinity of the exponential functions. Since the branch points come up in pairs, 
the leading term of the asymptotic expansion shows also a second branch point at x = 0. 
However, the behaviour of the function berx is regular at x = 0, therefore the asymptotic 
form in (P) is not suitable as an approximation for berx in the region near zero. In order to 
pursue the goal of picking up the right behaviour at infinity, but not to introduce undesired 



singularities in the zone of interest, a suitable auxiliary function has to be chosen. The 
auxiliary function should have the right ramification form at infinity, and the second branch 
point located outside the zone of interest (i.e., the negative axis) 

Only powers of the form An appear in the potential series of herx. In order to achieve 
the adequate efficency in the approximation, we should find such auxiliary functions and 
fractional approximations whose potential series only have exponents multiple of 4. 

Since the second branch point introduced with the asymptotic expansion must be outside 
the zone of interest, it is possible to choose as an auxiliary function Ai{x) = + x; 

however, that auxiliary function would not be efficient, since undesired power terms as 
x'^,x^, . . . would appear. A more suitable auxiliary function would be Ai{x) — + x^ 



or A2{x) — 1/v^l + T%^. However, in the last case the parameter r could be chosen in a 
convenient way in order to get better accuracy with the approximation, as we will discuss 
later, this can be considered as a free parameter. 

The other singularity that must be introduced through the auxiliary functions is of the 
form e^/^cos and e^^^sin The efficiency criteria limit the auxiliary functions 

by forcing them to have exponents of the form 4n. The above stated line of thought leads 
us to choose as convenient auxiliary functions 

cosh I I cos I I (4) 



and 

sinh (^) sin (^) 



x^ 



(5) 



Once the auxiliary functions are chosen as above, the fractional approximations auto- 
matically have only powers multiple of 4. 

All the previous considerations lead to the following form of the two-point quasifractional 
approximation for the real Kelvin function 

~ _ EUPkx'' cosh (^) COS (^) + ELo Pkx'' sinh (^) sin (^) 
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The auxiliary function vT+~a%^ has to be introduced in order to cancel the l/x^ be- 
haviour at infinity. With this function we also define a second free parameter a. 



III. CALCULATION OF THE TWO POINT QUASIFRACTIONAL 
APPROXIMATION TO THE REAL KELVIN FUNCTION 

Here we will consider only the simplest approximation to berx, thus n will be one and 
the approximation in eq. (P) will be reduced to 

(po + Pix") cosh (^) COS (^) + ^^i±f^(Po + Pxx^) sinh (^) sin (^) 

berx = ^y—L "^y^j. = \y—L ^y—!-. (7) 

(l + ga;4)^l + r%4 ^ ^ 

The ideas developed on the latest papers on quasifr actional approximations are in the sense of 
not to determine all the coefficients through the powers series and the asymptotic expansions 
but to leave one or two free parameters — a and r — that can be determined by minimising the 
maximum absolute error. It is clear that other methods to measure the discrepancy between 
the approximant and the exact function can also be used, such as the Lebesgue's integral 
of the quadratic difference, or other more elaborated methods using Lebesgue's integrals of 
the p-powers of the difference. However, for us, the most convenient way to measure that 
discrepancy has been the maximum absolute error. 

The main reason for using free parameters is to avoid the characteristic defects of classic 
Pade 10, also common in quasifractional approximations. Defects are when an extraneous 
pole appears near a zero in the numerator. That is, there is a zero in the numerator in xq 
and another zero in the denominator in x\ and the difference |a;o — X\\ is very small. Due to 
that problem, the approximation is in good agreement with the function in the whole zone 
of interest, except when close to xq and xi, where it goes to +00, —00 and 0. 

In order to determine the five unknowns po? -foiPi) -Pi; 9 three terms of the power series 
(I) and two terms of the asymptotic expansion (|^) will be used. So a linear system of five 
equations with five unknowns is obtained. The solutions will be given in terms of a and r. 

Looking at the exponential term in (^, the variable x/\/2 appears as more convenient. 



Thus the power series (0) will be written now as 

ber(xv^) = 1 - + (8) 

since we are using only three terms of this series we do not need to go further away than x^. 

Now, to compare the power series of ber x and ber x, we will first multiply both functions 
by (1 + gx^) in order to rationalise the right hand side. Furthermore the auxiliary functions 
appearing in (|^ and defined in @) and (|^), should be replaced by their respective power 
series. Thus we obtain 

(1 + -l.* + Jj^.») = 

(1 - irS' + ^^rV) ((1 - ix' + + p..') + 

where the first term in the first parenthesis on the right hand side is the power expansion 
of v^l + r*^x^, and constant factors have been absorbed in the definition of the parameters 
r and q. 

Equalising now the powers x* and x^ in both sides of the equation we obtain 

Po + ^^0 = 1 (10) 

bo + + -o?)Pq + -Pi + Vi - q = 

V 24 8 ^ 16 720 4 ^ " 2 ^ 16 



(11) 



'192 40320 12^ 



'5760 256 16 32 1440 3628800' 

(-- - -r')p^ + (--r^ - — + -a')P^ + -q=^ (12) 
^ 24 8 ^ 16 720 4 ^ 16^ 9216 ^ ' 

The leading term of the asymptotic expansion of ber x is 

[iPicosf^) +Piasin(^) 

berx~e^/^ ^ 

2 q^frx 

And comparing now the leading terms of ber x, we obtain 
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(13) 



A = ii^^ (15) 

V 27ra 



Equations (|T0]), (|TT]), (0), ([T^) and (p!5|) determine the values of po? -Po^Pi! -Pi aiid g as 
functions of a and r. 

When these equations are solved, an expression for each parameter is obtained in terms 
of a and r. And an approximation ber(x, a, r) is defined for each value of a and r. 

Given the initial values ao and tq, we numerically determine the value of the function 
|ber(x, a, r) — ber(x)| and we select the maximum value of this function, which will be the 
maximum absolute error. Therefore, the maximum absolute error, e(a;, r), will be a function 
of a and r. Looking now at e(a,r) as a function, we can determine the values am and 
which minimise this function as a two variable function. That procedure has been followed 
and the values a^. and are found as 

= 0.98 
Tm = 0.8367 

In order to avoid the defect in the approximation we have to consider only positive values 
of q, thus there is no problem with positive values for the variable x which is the region of 
interest. This must be taken into account when sweeping through a and r in the minimising 
procedure. Local minima which yield negative values for q must be discarded. 

Using these values of a and r in (p!OD, ([TT|), (|12]), (|1^ and (|15|), the parameters 



Po, Po,Pi, Pi and q are determined and the results are 

q = 27627.311660 
= -9750.649914 
Po = 19503.300340 
pi = 18628.544300 
Pi = 7873.669071 
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IV. ACCURACY OF THE TWO-POINT QUASIFRACTIONAL 
APPROXIMATION OF THE REAL KELVIN FUNCTION 

As it has been stated before, the discrepancy of the approximation will be determined 
by Aber x — ber x — ber x. 

A plot of the real Kelvin function, together with the approximation and ten times the 
absolute error is displayed in Figure 1. 

Besides, a table of the first five roots of the Kelvin function and the approximation is 
introduced. 

Notice that the highest relative error in the roots is the first one. Quasifractional ap- 
proximations to any function have normally their worst accuracy for values of order one. 
The behaviour of the approximation becomes more accurate as x increases. 

Another two measures of the good agreement of the approximation with the real Kelvin 
function is the relative difference of the position in x of the maxima of the function and the 
relative error of the amplitude of each maximum. The following table shows the position in 
X of the first five maxima of ber x and ber x as ber' x and ber x respectively, and the relative 
difference between them. The value of each of the first five maxima is also shown, appearing 
as ber^ x and ber^^ x for the real Kelvin function and the approximation respectively. The 
relative error of those figures also appears. 

Notice that all the relative errors, no matter what is being measured (i.e., the roots, the 
position of the maxima, the amplitude) decrease as x becomes greater. 

V. CALCULATION OF THE TWO-POINT QUASIFRACTIONAL 
APPROXIMATION FOR THE IMAGINARY KELVIN FUNCTION 

The same fine of thought carried out for the real Kelvin function will be followed in this 
section, since the form of the asymptotic expansions and the power series for both ber a; and 
heix are very similar. 
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The imaginary Kelvin function is defined as 0] 

herx = ^{J(){iVi x)) (16) 
where Jo{z) is the zeroth order Bessel function. The power series of heix is given by 

The asymptotic expansion of beix has been also given in several references. Again, our 
interest here lays mainly on the leading term, which can be written as 



e^/^(^sin(^)-VPeos(^) 

beix~ ^ ^= ^. (18) 

V27rx 

It is a characteristic of hypergeometric confluent functions to have an essential singularity 
at infinity. The asymptotic expansion picks that singularity up and shows it through a branch 
point at infinity with the exponential functions, in the same way that happened for ber x. 

The power series of beix goes as x'^,x^,x^^, .... In order to have an efficient approxi- 
mation, the auxiliary functions chosen for moving the artificially introduced branch point 
out of the positive axis and for introducing the singularities that appear in the asymptotic 
expansion should have a power series of the same form as heix. 

Following the same stream of thought used in determining the form of ber x, a correct 
form for the quasifractional approximation of bei x is 

^ ^ TT^n^oP^^'' cosh (^) cos (^) +SLo^.^^^ sinh (^) sin (^) 

where a and f are free parameters, introduced for the same reason explained for berx. We 
are going to use only the simplest form of the approximation, thus n = 1 will be substituted 
in the previous equation. Thus the two-point quasifractional approximation of bei x will 
give 

^(Pq+P^x^) cosh (^) cos (^) + (Po + Pix^) sinh (^) sin (^) 
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From now on, exactly the same procedure has been followed in determining the set of 
five equations used for finding the values of the parameters po, Po,pi, Pi, q in terms of a and 
f, and for minimising the discrepancy of the approximation. 

When dealing with quasifractional approximations it is common to take as an initial 
ansatz a and r equal to 1. This choice proved to be a valid starting point for berx, even 
though smaller differences are obtained when using the minimised Om and fm given in 
the previous section. The beix function is somewhat more pathological since the choice 
a = f = 1 yields a negative value of q, making the initial ansatz useless. 

For minimising the maximum absolute error of the two-point quasifractional approxi- 
mation of bei X a downhill simplex method has been used . The values obtained for the 
parameters are 





= 3.00 




= 3.00 


q 


= 19.11054940 


Po 


= -7.21235948 


Po 


= 15.42471896 


Pi 


= -30.32038957 


Pi 


= 24.39996523 



The plots of heix, heix and ten times the absolute error (lOAbeix = 10(beia; — heix) ) 
are shown in Figure 2. 

The relative error of the first root is smaller than the one of the first root of ber x because 
the value of the variable x is greater (i.e., the first root of beix it is not as close to zero as 
the one of berx). The accuracy for all the roots of berx and beix is very high in any case. 

Once again, the relative error related to the position of the first five maxima of the 
function and the error of their amplitude will be calculated in the same way performed for 
herx. The notation is exactly the same followed in the calculation of error for berx before 
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in Sec. 4, but adapted for the function beix. 



VI. ANALYTIC APPROXIMATION TO THE CURRENT DISTRIBUTION 
WITHIN A CYLINDRICAL CONDUCTOR 

The expression of the current distribution within a cyhndrical conductor due to a trav- 
elhng wave has been given in many references. The current distribution depends on the 
Kelvin real and imaginary functions and is given by 



J,(r) / hei\V2r/5) + bei'(v^r/5) \ 



J. 



(21) 



Vber2(v^ro/5) + bei2(v^ro/5) / 
where Jg is the current density on the surface of the conductor and tq is the radius on the 
surface. The skin depth is given by 5, and it depends on the frequency of the propagating 
wave. 

By means of substitution of berx and beix in (|2l|) it is possible to get an analytical 
approximation to the current distribution within a cylindrical conductor. 

The analytical approximation obtained has been used for studying the current distribu- 
tion within a copper wire with tq = 0.5mm. Four different values of the frequency have been 
taken. The same figures and plots could be found in Marion, for comparison purposes. 

The plots obtained using our approximation and those obtained using numerical com- 
putation are coincident and no difference can be found at this scale. Therefore, we are only 
showing the errors in Figure 3. 



VII. CONCLUSIONS 

In this paper analytic approximations have been found for the function bera; and heix. 
These have been used to calculate the current distribution within a cylindrical conductor 
produced by the propagation of one wave of frequency uj. 

In spite of the strong fluctuation characteristic of the functions berx and beix, the 
simple approximations found here not only reproduce the function with high accuracy for 
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all positive values of the variable but the infinite roots or zeroes of the function are also 
obtained and the relative error of each root decreases with the magnitude of the root. The 
largest relative error of the roots for bera; is about 2% and for beia; is 0.5%. 

The measure of the discrepancy is difficult for these functions because the function fiuc- 
tuates and, furthermore, the maximum of the amplitude increases and goes to infinity with 
X. Therefore, if we consider the absolute error this will go to infinity with x. 

On the other hand, if we study the relative error, these errors will go to infinity at the 
roots of the function. The best way to measure the good agreement of our approximation is 
just through the errors of the zeroes and of the amplitude of the oscillations. It is significant 
that these errors decrease with x despite of the large value of the function. The maximum 
error of the amplitude for berx is 3%, and 4% for heix in the value of the function and of 
0.28% and 0.80% in the site of the maximum for bero; and beia; respectively. 

These errors clearly show that our approximation can be used for most of the applica- 
tions we know until now, for example thermo-electrical corrosion in cylindrical conductors. 
In general, these approximations might be used whenever Kelvin functions appear and a 
numerical solution is not appropriate. As an example, we have applied our method to the 
simplest case of the current distribution in a cable, obtaining quite accurate results. 
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TABLES 



TABLE I. Comparison between the first five roots of berx and berx 


heix roots 






berx roots 


Relative 


error ( % ) 


2.84892 






2.78620 




2.2 


7.23883 






7.22030 




0.26 


11.67396 






11.66266 




0.094 


16.11356 






16.10548 




0.052 


20.55463 






20.54834 




0.032 


TABLE II. Comparison between the first five maxima of ber x 


and berx and their derivatives. 


her' X 


ber X 


err.(%) 


heVmX 


ber^x 


err.(%) 


6.03871 


6.0215 


0.28 


-8.86404 


-8.61484 


2.9 


10.51364 


10.5027 


0.10 


153.782 


151.237 


1.7 


14.96844 


14.9605 


0.053 


-2968.68 


-2933.93 


1.2 


19.41758 


19.4114 


0.032 


60161.2 


59616.8 


0.92 


23.86430 


23.8592 


0.021 


-1.25374 X 10^ 


-1.2445 X 10^ 


0.75 


TABLE III. Comparison between the first five roots of heix and heix. 


heix roots 






heix roots 


Relative 


error ( % ) 


5.02622 






4.99873 




0.55 


9.45541 






9.44110 




0.15 


13.89349 






13.88400 




0.069 


18.33398 






18.32689 




0.037 


22.77544 






22.76977 




0.024 
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TABLE IV. Comparison between the first five maxima of heix and beix and their derivatives. 



bei X 


bei X 


err.(%) oeimX 


beima^ 


err.(%) 


3.77320 


3.74307 


0.80 2.34615 


2.259 


3.9 


o.zouyy 


o.zoi io 


U.Ll -OD.iD04 


-oo.4Uoy 




12.74215 


12.7329 


0.073 670.16 


660.955 


1.4 


17.19343 


17.1865 


0.040 -13305.5 


-13169.6 


1.1 


21.64114 


21.6356 


0.026 273888 


271661 


0.82 


TABLE V. Skin depth at different frequencies. 


Case 




u = hjjlv: 




b (mm) 


1 




10^ 




2.1 


2 




10^ 




0.66 


3 




10^ 




0.21 


4 




10^ 




0.066 
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FIGURES 
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FIG. 1. Kelvin real function and its two-point quasifractional approximation, along with ten 
times the discrepancy between them. 
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FIG. 2. Kelvin imaginary function and its two-point quasifractional approximation, along with 
ten times the discrepancy between them. 
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FIG. 3. Plots of the discrepancy between the analytical current distribution within a cylindrical 

conductor and the results obtained with two-point quasifractional approximation for four different 
values of the skin depth. 
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